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Abstract. The auxiliary field method is a technique to obtain approximate 
closed formulae for the solutions of both nonrelativistic and semirelativistic 
eigenequations in quantum mechanics. For a many-body Hamiltonian describing 
identical particles, it is shown that the approximate eigenvalues can be written 
as the sum of the kinetic operator evaluated at a mean momentum po and of 
the potential energy computed at a mean distance rg. The quantities po E'nd tq 
are linked by a simple relation depending on the quantum numbers of the state 
considered and are determined by an equation which is linked to the generalized 
virial theorem. The (anti)variational character of the method is discussed, as well 
as its connection with the perturbation theory. For a nonrelativistic kinematics, 
general results are obtained for the structure of critical coupling constants for 
potentials with a finite number of bound states. 
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1. Introduction 

The auxiliary field method (AFM) is a very powerful method to obtain approximate 
analytical expressions for the eigenvalues of one, two and many-body systems with 
both nonrelativistic or semirelativistic kinematics. It has been shown in a series 
of paper [TJ[2l[3l|4j[5l|6l[7l|8] that it can be applied with great success in many 
physical situations. The basic idea is to replace a problem which is not solvable, 
for example because of a complicated potential or a semirelativistic kinematics, by 
another one which can be treated analytically. In so doing, it is necessary to introduce 
auxiliary fields Uf^. The original Hamiltonian H is replaced by a new Hamiltonian 
H(pk), called the AFM Hamiltonian. If these auxiliary fields are chosen as i>fc(0) 
in order to extremize the AFM Hamiltonian, this one coincides with the original 
Hamiltonian: H{i>k(0)) ~ H. Thus, both formulations are completely equivalent. 
The approximation lies in the fact that the auxiliary fields are considered no longer as 
operators, but as a real constants Vk- An approximate value of the exact eigenenergy 
E is then given by an extremal eigenenergy E{i'k{0)) of the AFM Hamiltonian 
H{i'k), which is in principle much simpler than H. An approximate state for the 
corresponding eigenvalue can also be obtained. The quality of this approximation 
has been studied and discussed in detail in the papers mentioned above. Among the 
interesting properties of the AFM, we can mention its great simplicity and its ability to 
treat on equal footing the ground state and the various excited states. This procedure 
was first introduced to get rid off the square root kinetic operator in calculations for 
semirelativistic eigenvalue equations [9j [10] . As the AFM is an extension of these first 
calculations, we just keep the same name for the method. 

As it is shown in [TT] , the AFM has strong connections with the envelope theory 
[121 [131 [m [H] ■ Nevertheless both methods have been introduced from completely 
different starting points. In particular, the AFM introduces the notion of auxiliary 
fields which is a key ingredient to interpret the method as a mean field approximation 
and which can be very useful to compute mean values of observables [1] [^ . 

In this work, we present some new and general properties of the AFM. So, only 
the basic ingredients necessary for the understanding of the subject treated here are 
recalled in Sect. [31 We refer the reader to our works mentioned above or to our 
review paper [16 for an exhaustive overview of the method and its applications. New 
results or generalizations are presented in the following sections. The connection of 
the AFM with the generalized virial theorem is presented in Sect. [3l The use of 
the perturbation theory for the AFM is explained in Sect. [H For a nonrelativistic 
kinematics, the general structures of critical coupling constants for potentials with a 
finite number of bound states are presented in Sect. [SI A summary of the results is 
given in Sect. [H 

2. The auxiliary field method 

Let us consider a system composed of N particles, interacting via one-body potentials 
Ui and two-body potentials Vij, and moving with a nonrelativistic or a semirelativistic 
kinetic energy. In principle, the AFM can treat this problem in such a general form, 
but it is manageable in practice only if the particles are identical. This implies that 
they all have the same mass to, that the form of the one-body potentials is the same 
for all particles Ui = U, and that the form of the two-body potentials is the same for 
all pairs of particles Vij = V . So, the most general Hamiltonian we will consider in 
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this paper has the following form 

N N N 

H = J2^P^ + m^ + Y.Ui\s,\)+ J2 ^(l^^.l). (1) 

i—l i—1 i<j — l 

with Si — Ti — R and rij = ri — Tj. is the position of the particle i, Pi is its 
conjugate momentum and R is the position of the center of mass of the N particles 
[N > 2). It is assumed that X^iLi — 0. Following the AFM, this Hamiltonian is 
"replaced" by the auxiliary Hamiltonian iJ, with auxiliary potentials P{x) and S{x), 
and which depends on 3 auxiliary fields /i, v and p 

i=l " i=l i<j = l 

provided the states considered are completely symmetrized. Most of the results 
presented in this section come from [B] , but we use here the more convenient notations 
developed in [SI US]- The function B{^,i/,p) is not useful to detail in this work, but 
it can be rebuilt from results in [B]. If this Hamiltonian is analytically solvable, an 
AFM analytical approximation of a mass M of the A^-body Hamiltonian H is given 
by an eigenvalue AIq of the Hamiltonian H(po,vo, po) where the 3 optimal auxiliary 
parameters po, and po extremize this eigenvalue. These parameters depend on the 
quantum numbers of the state. 

At this stage, it is not obvious that the solution Mq is a good one. But, the 
comparison theorem of the quantum mechanics can be used to obtain significant 
information about the AFM eigenvalues. This theorem states that, for some eigenvalue 
equations, if two Hamiltonians are ordered, T?^^-* < T?*-^-* ((T?^^-*) < (iJ^^^) for any 
state), then each corresponding pair of eigenvalues is ordered -Ejgj < -^'{ej.i where 
{0} represents a set of quantum numbers. This inequality can be obtained from the 
Ritz variational principle [17 , but it can also be derived from the Hellmann-Feynman 
theorem [18) . If we can show that the auxiliary Hamiltonian H{pq,vq^ p^) is greater 
or lower than the genuine Hamiltonian iJ, then it is possible to use the comparison 
theorem to locate the AFM eigenvalues with respect to the exact ones. 

In the case of a nonrelativistic kinematics. The AFM yield an upper (lower) bound 
if the potentials U {x) and V{x) could be bounded from above (below) by the auxiliary 
potentials P{x) and S{x) respectively [6]. For a semirelativistic kinematics, the AFM 
implies a replacement of the square root operators by a nonrelativistic form of the 
kinetic energy (see ^) and this yields to an increase of the eigenvalues [5J[TS]. So, 
in this case, the AFM solutions are upper bounds of the exact ones if the potentials 
U {x) and V{x) can be bounded from above by the auxiliary potentials P{x) and 
S{x) respectively. In other cases, nothing can be said about the possible variational 
character of the solutions. 

Let us note that a lower bound for the ground state (and then for the whole 
spectrum) of the general Hamiltonian ^ for a boson- like system has been proposed 

X A boson-like system is composed of particles whose total spatial wavefunction can be completely 
symmetrical. For instance, this is the case for a system of quarks inside a baryon: quarks are fermions, 
but the baryon is characterized by a completely antisymmetrical colour function so that the rest of 
the total wavefunction must be completely symmetrical. Similarly, a fermion-like system is composed 
of particles whose total spatial wavefunction can be completely antisymmetrical. 
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in [IH] . It takes the following general form (the one-body potential is introduced using 
the fact that |si| = jsal = |ri2|/2 for iV = 2) 



M > N inf 



U 



1, 



N-1 



V{\r\) 



(3) 



but also works for nonrelativistic kinematics (20j. In this latter case, using the AFM 
results with = 2, a lower bound for the mean value in ^ can be computed provided 
the potentials U{x) and V{x) could be bounded from below by the auxiliary potentials 
(see Sect.[2J]). 



2.1. The case N >2 

For arbitrary values of N, the Hamiltonian ([5]) is entirely analytically solvable for 
the unique choice P{x) and/or S{x) equal to x'^. It can then be shown [5] that the 
non linear system determining the 3 variables (/io, Po) can be recast in the form of 
one transcendental equation depending on the single variable ~ ^y2f^o{'^o + Npo). 
Moreover, an eigenmass can be computed from the Xq quantity only. Thus, the 
eigenvalue problem for the TV-body system can be determined simply by the set of the 
following two equations: [5] 




(4) 



(5) 



where Q is a global quantum number (see below), where K{x) = U'{x)/P'{x) = 
U' (x) / {2x) and L{x) — V (x) / S' (x) = V'{x)/{2x), and where the number of pairs 

N{N- 1) 



Cn 



(6) 



has been introduced for convenience. The prime denotes the derivative with respect 
to the argument. In this framework, an approximate AFM eigenstate is given by an 
eigenstate of H{fj,Q,VQ, pq). It is written in terms of Jacobi coordinates as a product 
of (A^ — 1) oscillator states with sizes depending on A^ and Xq A nonrelativistic 
version of (01)- © can be obtained in the limit m — oo In this case, po — >■ m. 
Further simplifications occur also for the ultrarclativistic limit m = 0. 

A state depends on (A^ — 1) radial quantum numbers nt and (A^ ~ 1) orbital 
quantum numbers h, as well as intermediate coupling quantum numbers which are 
not considered here. The global quantum number resulting from the AFM treatment 
is then 



N-l 



(7) 



All quantum numbers are not allowed, depending on the nature of the particles. In 
particular, the ground state for a boson-like system is just Q = 3{N — l)/2, while the 
ground state of a fermion-like system is much more involved and needs the introduction 
of the Fermi level [Q. 

Since P{x) and/or S{x) equal to x^ can only be used, an upper bound is computed 
for most of the relevant interactions, a fortiori for a semirelativistic kinematics. For 
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instance, an AFM mass formula has been obtained for a system of N relativistic 
massless quarks interacting via a hnear one-body confinement and a two-body 
Coulomb potential (this kind of Hamiltonian is pertinent for variant theories of the 
quantum chromodynamics). The accuracy of this formula has been numerically tested 
in [6] with TV = 3: Relative errors less than 20% have been obtained for the lowest 
states. It has also be shown in [21] that the A''-dependence of this formula is the correct 
one for N ^ oo. When a closed formula cannot be computed, numerical solutions 
(generally upper bounds) can always be easily obtained. This is valuable for a iV-body 
system. 

2.2. The cases N and N ^ 2 

For iV = 2, Si = — S2 = ri2/2. So, the potential U{x) becomes redundant with 
the potential V{x) and can be ignored. Moreover, the Hamiltonian H simplifies 
because pi = ~P2 = P- Thus, for a nonrelativistic kinematics, the case of two 
different particles can be considered by replacing the kinetic part 2m + p^/m by 
mi -|- m2 +p'^/{2mr) where rrir is the reduced mass. A priori, above calculations are 
only valid for N > 2. But, starting from the one-body equivalent of Hamiltonian ([T]), 
it can be shown that equations (31)-® are also relevant for = 1 by setting V{x) — 
and reinterpreting pi and Si as conjugate variables. 

For both A^ = 1 and 2 systems, the more general form sgn(A) x'^ can be used for 
the auxiliary potential, instead of only x'^. This leads to various expressions for Q. The 
complete calculation shows that the same system (II])-® is found and that the only 
trace of the auxiliary potential lies in the structure of the global quantum number Q. 
In practice, Q = 2n + /+3/2 with P{x) or S{x) = x^ (see dTj with A^ = 2), Q = n+l + 1 
with P{x) or S{x) = -l/x [3] and Q = 2{~an/3f^^ for S-wave states with P{x) or 
S{x) = X [7], where an is the (n -I- 1)^'^ zero of the Airy function Ai. Depending on 
the kinematics, closed form formulae have been obtained for various potentials: sum 
of two power-law, logarithmic, Yukawa, exponential, square- root [U [21 [3l SJ [3 [6] . If 
a closed formula cannot be computed, the method is then not really interesting since 
a lot of numerical techniques can be harnessed to find accurate solutions for one- or 
two-body systems. 

For some nonrelativistic systems, it is possible to use two forms of the auxiliary 
potential to obtain both upper and lower analytical bounds of the exact solutions. 
The following potentials, aln(6a;), ax — b/x, or sgn(A) ax^ (with a > 0, 

6>0, — 1<A<2) can be bounded from below (above) with the choice —1/x (x^) 
for the auxiliary potential. For instance, let us consider the Hamiltonian 

H=^ + Va\^+P. (8) 
The eigenenergies computed with the AFM gives 

E^,u = ^{ouy) + ^) With y-'j{§^f. (9) 

and where G-{Y) is the solution of the equation 4G_(r)'* ~ 8G-{Y) - 3Y = 0. 
Upper (lower) bounds are obtained with Q — 2n + l + 3/2 [Q — n + l + \). The quality 
of these bounds are studied in [H [5] where more details are given about this solution 
and the function G-{Y). 



Further developments for the auxiliary field method 



6 



3. Connection with the virial theorem 

The general virial theorem links the mean values of the directional derivatives of the 
kinetics operator and the potential [22l [23] . Using the Hellmann-Feynman theorem as 
in [23], it can be applied to the general TV-body Hamiltonian ([T]) to yield 

N {pk ■ Vp,T{pk)) =N{si- Vs,U{si)) + Cn {n, ■ V,,^F(n,)) , (10) 

with arbitrary numbers {fc, 1^ i ^ j} if the mean values are taken with a completely 
symmetrized eigenstate of the A^-body Hamiltonian. The operator T is defined by 
T{x) = y/x^^+m? or by its nonrelativistic counterpart m + x^ / {2m). Let us introduce 
the distance tq = NQ / Xq and the momentum Pq — Q/ro- It is a simple algebra 
exercise to show that formulae (j!])-® can be written as: 

Mo^NT{po)+Nu[^)+CNv(^-^y (11) 
Po = (12) 

,13, 

These equations have not been presented in our previous papers. Before discussing 
their physical meaning, let us look at the quantities rp and pq. Using formulae of the 
appendixes in I?] and [2T], the following observables can be analytically computed: 

N \ 

(14) 




(15) 

This shows that ro can be considered as a mean radius for the system and po as a 
mean momentum per particle. Indeed, (jl4p and (|15p imply that 



^)=Po, v^=^, y^=^, (16) 

for arbitrary i ^ j since the mean values are taken with completely symmetrized 
states. These results can also be obtained using the more general relations (66)-(68) 
in [B] relevant for P{x) and S{x) different from x^. 

With this new formulation, an AFM eigenvalue given by (jlip is simply the kinetic 
operator evaluated at the mean momentum po plus the potential energy computed at 
some mean radius depending on tq. As one could expect, the kinetic energy and the 
one-body potential energy are proportional to the number of particles and the two- 
body potential energy is proportional to the number of pairs. Formula ()lip looks like 
a semiclassical approximation but this is absolutely not the case. The AFM yields an 
approximate A^-body wavefunction [6l [21] , and the relation (fT2|) between pa and ro 
is a full quantum link, function of the quantum numbers of the system. At last, the 
value of ro (and thus of po) is the solution of a transcendental equation (fT3]) which is 
the translation into the AFM variables of the generalized virial theorem (flUl) which 
comes from very general properties of quantum mechanics. These considerations prove 
that the AFM really relies on very sound physical basis. Once the system (|lip - p^ is 
written, it can appear finally quite natural to obtain such a result. The problem is to 
find a relevant link between the mean values ro and po- This is solved by the AFM. 
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It is generally possible to improve the quality of the AFM eigenvalues with a 
slight modification of the principal quantum number. A particularly simple form 
which works quite well is given by 

Q= ^(an, +7(^-1), (17) 
1=1 

where the values of parameters a, /3 and 7 depend on both the interaction and the 
kinematics. They can be determined by an analytical procedure in some cases by 
using analytical results coming from WKB approximations or variational calculations 
P,, 8 . Even if it less interesting, a fit on numerically computed exact eigenvalues can 
always be implemented [T] [5J [31 131 [3] . With the form p7)) . the variational character 
of the AFM approximation is lost, but the relative errors can be sometimes strongly 
reduced. 



4. Connection with the perturbation theory 



It has been shown in l^, that, for one- and two-body nonrelativistic systems, the 
AFM and the perturbation theory give similar results when the potential is an exactly 
solvable one plus a small perturbation. This result is extended here for the general 
Hamiltonian ([IJ, that is to say: N particles, semirelativistic kinematics and arbitrary 
potentials U{x) and V{x). 

Let us first assume that each pairwise potential V^d^'ijl) is supplemented by a 
term e?;(|rij|), with e ^ 1 in order that ev{x) <Si V{x) in the physical domain of 
interest. In the system pT|) - (IT3)) . the potential V{x) is replaced by V{x) + ev{x). 
In this case, new values ri and pi for the mean radius and momentum will be the 
solution of the new system: 



Ml 



Pi n 



NT{pi)+Nu(^^^+Cn 



V 



ri 



'Cn 



Np,T'{p,) = r^U' (^^ 



V 



ICn 



ri 



(18) 
(19) 
(20) 



In this case, power 
(1 - 6)pq from 
T'ipo) - SpoT"{po), 



'Cn , 

Writing ri = (1 -I- (5)ro, we can expect 5^1 since e ^ 1. 
expansions at first order can be computed. We have p\ ~ 
and we can write T(pi) ss T{po) — SpoT'{po), T'{pi) 

U{ri/N) « U{ro/N) + 6ro U'{ro/N)/N, etc. Equation ^ reduces to an expression 
of the form 5 « e /i(ro) where /i is a quite complicated function of T', [/', V' and their 
derivatives. It confirms that 6 ^ 0(e). The precise form of h is given below in the 
most general case. It is then possible to perform an expansion of Mi to obtain 

'ro- 

^0 \ , /TT- , -,r! I \ ^ ^ f ro 



N5poT'{po) 

C^droV 



'Cn J \ \JCn 

Using pT|) and (|13p . this equations simplifies to 




Nu(^-^)+5roU' 

+ Cn (-V 



'Cn 



+ 0(6^). 



Ml = Mo + Cn€v 



(21) 



(22) 



This result could seem quite obvious, but it demonstrates that the knowledge of tq is 
sufficient to obtain the contribution of the perturbation at the first order. 
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Let us now look at the most general case and assume too that each [one-body 
potential C/(|sj|) / kinetic operator r(|pj|)] is supplemented by a term [77u(|si|) / 
Tt{\pi\)], with [77 < 1 / r < 1] in order that [r]u{x) < U{x) / Tt{x) < T{x)] in the 
physical domain of interest. With similar calculations, we finally find 

Ah = Mo + NTtipo) + N7^u(^^^+CNev (^^^^ + 0{e^ ,rj^T^). (23) 

The parameter 5 is determined at the same order by the following relation 

NpQTt' (po) -rovu' (^^^ - 



Cn 

2Np,T' (po) + NplT" (po) + |c/" (^) + rl V" 



(24) 



Perturbed observables and wavefunctions can then be computed at first order, since 
ri = (1 + (5)ro and pi = (1 — 5)pq at this order. 

The contribution of a perturbation at the first order can thus be very easily 
computed within the AFM once the unperturbed problem is solved. In order 
to check the quality of this approximation, let us consider a case in which the 
unperturbed Hamiltonian H can be solved exactly by the AFM, that is A/q is the 
exact solution. If the small perturbation potential is written e X]i<j=i ^(l^y Dj 
quantum perturbation theory says that the solution is given by 

= A/o + Cn t {v{\r,,\)) + 0(6^), (25) 

for any pair (ij). The mean value is taken with a completely symmetrized eigenstate 
of the unperturbed Hamiltonian H. The comparison of (1^51) with ([22) shows that 
(u(|rij|)) is replaced by v [r^j \JC-n) within the AFM. This is to be compared 
with the exact relation (5(|ry|)) = S [r^j \JC-n) for the auxiliary potential [HIS]. 
So, the AFM does not give the same result as the perturbation theory. But the 
agreement can be very good, as shown with several examples calculated explicitly in 
[2]. Similar discussions can be made for small one-body perturbation potentials or 
small perturbations of the kinematics. 



5. Critical coupling constants 

Some interactions, as the Yukawa or the exponential potentials, admit only a finite 
number of bound states. Let us assume that such an interaction can be written as 
Wix) — —K,w{x), where k is a positive quantity which has the dimension of an energy 
and w{x) a "globally positive" dimensionless function such that lima;_^(x) w{x) = 0. 
We can introduce the notion of critical coupling constant k{{6}) where {6} stands for 
a set of quantum numbers. This quantity is such that, if k > k{{9}), the potential 
admits a bound state with the quantum numbers {9}. The interaction energy for 
the state with quantum numbers {9} is then just vanishing for k — k({9}). We refer 
the reader to [351 UHl 1271 HE] for detailed explanations about how to compute critical 
coupling constant in a given potential. 

Let us consider a nonrelativistic A^-body system (no manageable calculation can 
be performed for a semirelativistic kinematics) with one-body potentials U{x) = 
—ku{x) and two-body potentials V{x) = —gv{x), both independent of the particle 
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mass and both admitting only a finite number of bound states. Tlie system pT |) -p ^ 
for a vanisliing energy gives: 



^9 2 

2 TO 



NkNuQ)+CN9Nvi^-^j , (26) 

- kNTou' (^^^ - ^/C^ gNrov' (^-^^^ , (27) 

where and are the critical constants for the system with TV particles. The 
elimination of the ratio N Q'^/{mr'^) from both equations yields the equality 

2N kN u (^) + 2Cn gN V (^-^ 

= -fcArrou' (^^^ - yC^gArro-u' ^-^^^ . (28) 

When potentials u and v are both taken into account, nothing interesting can be said. 
So let us consider one type of potential at once. 

Assuming that only two-body forces are present, ([SS)) reduces to 

where the parameter g^ has disappeared. Introducing the new variable yo = fo/VCN, 
we can rewrite (PSI) and (l?f)) as: 

1 2 

^^^yg^(yo) NiN-l)^m' ^^^^ 
2viyo)+yQv'iyo) = 0. (31) 

The variable yo, determined by (j3ip . is independent oi N , Q and m, and depends only 
on the form of the function v{x). So, the general formula (jSOp . which was not obtained 
in our previous works, gives precise information about the dependence of the many- 
body critical coupling constant ^at as a function of all the characteristics of the system. 
With the system (p0| -(|3i p . it is easy to recover some limited previous AFM results 
obtained for the critical coupling constants of Yukawa and exponential interactions 
[21 16]. For instance, with the two-body Yukawa interaction V{x) = —g exp{~j3x)/x, 
we have 

2e/3p2 
N{N- Xfm' 

For N = 2 and Q = n -I- / -t- 1, reasonable upper bounds of the exact critical coupling 
constants are obtained f^. 

Within the AFM approximation, the ground state (GS) of a boson-like system 
is characterized by Q = |(iV — 1). We obtain in this case the following very general 
relation valid, at the AFM approximation, for all pairwise potentials with a finite 
number of bound states 

gjv+i(GS) ^ 

5iv(GS) N + \ ^ > 

This ratio has previously been obtained and numerically checked for several 
exponential- type potentials [291 130| . Similarly, in the same general situation, 

5^(GS) = ^52(GS), (34) 
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indicating that in order to bind a A^-body system, a coupling N/2 times smaller than 
the coupling for a two-body problem is sufficient [HI [30] ■ 

Assuming that only one-body forces are present, a similar calculation gives: 

2u{yo) + you'{yo)^0, (36) 

where the change of variable yo = ro/N has been used. Again, the general formula 
(|35p . which was not obtained in our previous works, gives precise information about 
the dependence of the one-body critical coupling constant /cjv as a function of all the 
characteristics of the system. For the ground state of a boson-like system, we obtain: 



fcw+i(GS) / 



2 

(37) 



fcA'(GS) =a(^^^^ fc2(GS). (38) 

These results are strongly different from those for pairwise forces. 

If the AFM gives upper (lower) bounds for the exact eigenvalues, the critical 
coupling constants predicted by formulae above are upper (lower) bounds for the 
exact critical coupling constants. 



6. Summary 

The main interest of the auxiliary field method is to obtain approximate closed 
formulae for the solutions of nonrelativistic and semirelativistic eigenequations in 
quantum mechanics. The idea, strongly connected with the envelope theory, is to 
replace a Hamiltonian H for which analytical solutions are not known by another one 
H which is solvable and which includes one or more auxiliary real parameters. The 
approximant solutions for H, eigenvalues and eigenfunctions, are then obtained by the 
solutions of H in which the auxiliary parameters are eliminated by an extremization 
procedure for the cigcnenergies. The AFM can yield upper or lower bounds (both 
in some favorable situations) on the exact eigenvalues. The nature of the bound 
depends on the fact that H > H or H < H. With a semirelativistic kinematics, only 
upper bounds can be obtained because of the replacement of the kinetic operator by 
a nonrelativistic one. For many-body systems, only one type of Hamiltonian H can 
be used. So, it is not possible to obtain both upper and lower bounds for the whole 
spectrum in this case. Nevertheless, for a nonrelativistic kinematics, a lower bound 
for the ground state can be sometimes computed. 

Provided the structure of the Hamiltonian H is well chosen (nonrelativistic 
kinematics plus power- law potentials) , an eigenvalue computed by the AFM is simply 
the kinetic operator evaluated at a mean momentum pq plus the potential energy 
computed at some functions of the mean radius rp. The product rgpo is equal to a 
global quantum number characterizing the state considered, and the value of rg (and 
then oi po) is the solution of a transcendental equation which is the translation into 
the AFM variables of the generalized virial theorem. This new result gives sound 
physical basis to the method. 

Once a problem is solved within the AFM, it is very easy to compute the 
contribution of a small perturbation at the first order. It is given by the perturbation 
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Hamiltonian evaluated at the mean momentum pq for a kinetic energy or at a function 
of the mean radius tq for a potential. The resiilt does not coincide with the one 
obtained by the quantum perturbation theory, but the agreement can be very good. 

The AFM gives a very general formula for the critical coupling constants of 
nonrclativistic Hamiltonians with a finite number of bound states. The dependence 
on the quantum numbers, the mass m of the particles, the number N of particles, 
and the structure of the potential are predicted. Different N behaviours are obtained 
depending on the one-body or pairwise character of the interaction. If the AFM 
gives upper (lower) bounds for the exact eigenvalues, the critical coupling constants 
predicted are upper (lower) bounds for the exact critical coupling constants. 
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